Strain gradient induced polarization in SrTiOa single crystals 
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Piezoelectricity is inherent only in noncentrosymmetric materials, but a piezoelectric response 
can also be obtained in centrosymmetric crystals if subjected to inhomogeneous deformation. This 
phenomenon, known as flexoelectricity, affects the functional properties of insulators, particularly 
thin films of high permittivity materials. We have measured strain-gradient-induced polarization 
in single crystals of paraelectric SrTiOa as a function of temperature and orientation down to and 
below the 105 K phase transition. Estimates were obtained for all the components of the flexoelectric 
tensor, and calculations based on these indicate that local polarization around defects in SrTiOa may 
exceed the largest ferroelectric polarizations. A sign reversal of the flexoelectric response detected 
below the phase transition suggests that the ferroelastic domain walls of SrTiOs may be polar. 



Flexoelectricity is the coupling between dielectric po- 
larization and strain gradient. Because strain gradients 
break inversion symmetry, flexoelectricity allows extract- 
ing charge from deformations even in materials that are 
not piezoelectric. The historical development of this phe- 
nomenon has been summarized in the reviews by Tagant- 
sev [l| and more recently by Cross [3] and Sharma et al. 

For decades after its proposal in the early 1960s [HQ, 
there were very few studies of flexoelectricity, which may 
seem surprising given that it is a more general property 
than piezoelectricity. The reason is the comparatively 
small magnitude of the flexoelectric coefficient / ~ e/a 
(where / is the constant of proportionality between po- 
larization and strain gradient, e the electronic charge and 
a the lattice parameter) in typical dielectrics which, com- 
bined with the small size of strain gradients in bulk sam- 
ples, means that the flexoelectric effect is generally mi- 
nuscule. Two recent developments, however, are chang- 
ing the perceived importance of flexoelectricity. First, 
a predicted proportionality between flexoelectric coeffi- 
cient and dielectric constant has inspired Cross and 
coworkers to conduct a series of seminal experiments in 
which they have measured the flexoelectric coefficients of 
high permittivity materials such as relaxors and ferro- 
electrics [3 0] , and found them to be very large indeed. 
In addition, Fousek et al. 8] realized that this large 
flexoelectric response could be used to produce piezo- 
electric composites from centrosymmetric materials by 
simply tailoring their shapes, an idea which was recently 
demostrated experimentally by Zhu et al. Second, 
recent works by Catalan et al. [Io| have shown that fer- 
roelectric thin films can sustain large strain gradients due 
to lattice mismatch with the substrate, and that the flex- 
oelectric effect of such gradients can have an important 
effect on their functional properties, particularly the fre- 
quently discussed lowering of permittivity in ferroelectric 
thin films. 

In spite of its growing relevance, our understanding of 
flexoelectricity is still rather incomplete. On a theoret- 
ical level, there are still no first-principles studies. On 



an experimental level, only ceramic samples have been 
analyzed, with the problem that grain boundaries can 
contribute to measured charge due to their possibly po- 
lar nature [III, through tribological effects or via sur- 
face piezoelectricity Also, because of the need for a 
high dielectric constant, only materials with ferroelectric 
phases (or polar nanoregions in the case of relaxors) have 
been studied so far. This has the drawback for interpreta- 
tion that ferroelectricity (and thus piezoelectricity) may 
persist above the nominal phase transition temperature 
due to local strain effects, and contribute piezoelectrically 
to the measured charge. 

We have chosen to work instead with single crystals 
of the paraelectric strontium titanate. At room tem- 
perature SrTiOa (STO) has the cubic perovskite struc- 
ture and remains centrosymmetric even in the tetragonal 
phase below 105 K. The absence of piezo- and ferroelec- 
tric contributions, combined with relatively high dielec- 
tric permittivity, makes STO a natural choice for study- 
ing flexoelectricity. The use of single crystals of different 
orientations also allows all the flexoelectric tensor com- 
ponents of STO to be determined. 

High purity STO single crystals of various dimensions 
(typically 3-5mm wide, 5-15mm long and 50-500/im 
thick) were obtained from PiKem Ltd. Impurity levels (in 
parts per million) determined by the supplier are: Ni=3, 
Fe=2, Cr<2, Ba, Na and Si each <1. Ca impurities at 
<lppm are several orders of magnitude below the criti- 
cal concentration of Ca 2+ ions above which a ferroelec- 
tric phase transition can be induced at low temperatures 
[13]. The transparent colourless crystals were supplied 
with a surface roughness of a few A. Top and bottom 
Au electrodes of area 10-30mm 2 (depending on crystal 
size) were deposited by sputtering. Pt wires (50/xm in 
diameter) were attached to the Au electrodes with sil- 
ver paste which was then annealed at 130°C to improve 
conductivity and mechanical robustness. 

The experimental setup is sketched in figure [T] A dy- 
namical mechanical analyzer (DMA) with an insulating 
quartz probe was used to induce oscillatory bending (typ- 
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ically driven at 30-40Hz) and measure its amplitude. A 
static stress was applied simultaneously to hold the sam- 
ple in place. For temperature dependent measurements, 
heating and cooling could be achieved by competitive 
action of a resistive heater and a liquid N2 bath. The 
displacement current / due to the induced polarization 
was measured using a Signal Recovery 7265 dual phase 
lock-in amplifier. The maximum strains (approx. f CP 6 ) 
achieved by the bending during the experiments were far 
below those which, according to thermodynamic [3, [3| 
and ab-initio [ll| calculations, are capable of inducing a 
polar phase in the investigated temperature range. The 
crystals are therefore expected to be neither ferro- nor 
piezoelectric. A direct measurement of piezoelectricity 
showed no signal above the noise level placing an upper 
limit on the piezoelectric coefficient ^33 of 0.03pC/N. In 
addition, no signs of second harmonic generation (SHG) 
could be detected. 
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FIG. 1: Experimental setup for flexoelectric measurements. 

The strain gradient was derived from the usual equa- 
tion for a bent beam (e.g. see as 
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where L is the distance between the knife edges (in our 
case L =10, 7.5, or 5mm), zq is the displacement at the 
centre as measured by the DMA and the distances Xi are 
measured from the centre of the crystal. 

Flexoelectricity is described by a fourth rank tensor 

fijkl 
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where P is the flexoelectric polarization and eki is the 
symmetrized elastic strain tensor. If tu/2ir is the fre- 
quency of the applied mechanical stress and A is the 
electrode area, the average out of plane polarization can 
be computed by measuring the ac current produced by 
bending and using P3 = I/wA. The effective flexoelectric 
coefficient (i was calculated from the measured average 
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FIG. 2: Temperature dependence of the ratio I/zo (which 
is proportional to fi) and of the capacitance for a sample 
polished down to 100/im using diamond impregnated lapping 
film. The flexoelectric response increases with the dielectric 
constant, showing an anomaly at the ferroelastic phase tran- 
sition. The absolute temperature values shown may be a few 
degrees too low due to some thermal lag between the thermo- 
couple and the sample. Inset shows the linear dependence of 
flexoelectric polarization on strain gradient for 300/im thick 
STO crystals of different orientations at room temperature. 



polarization and strain gradient using 
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where a is the half-length of the electrodes. 

From phenomenological arguments, the flexoelectric 
coefficient is expected to be proportional to the dielectric 
constant To test this, the flexoelectric current was 

measured as a function of temperature (figure [5]). As ex- 
pected, the current increases upon cooling, qualitatively 
following the dielectric permittivity. On approaching the 
ferroelastic transition, however, an anomaly, not present 
in the dielectric constant, is seen in the flexoelectric re- 
sponse. 

Previous studies 



17, UM have shown that ferroelas- 



tic domains appear below 105K and are responsible for 
the large softening of the elastic modulus in STO. The 
relaxation of strain gradients due to domain wall read- 
justments may therefore be expected to reduce the flex- 
oelectric response, as indeed observed. This gradient 
relaxation by domain readjustment can be studied by 
changing the static force applied to the crystal. With 
large static force, the crystal is very bent and the do- 
main walls impinge on each other, so that they can no 
longer move under the dynamic load. Thus, the mechan- 
ical and flexoelectric response should approach that of a 
monodomain crystal. For low static forces, the domain 
walls are quite free to move and thus the relaxation of 
the strain gradient is maximum. Accordingly, one would 
expect the flexoelectric current to drop. Instead, how- 
ever, we observed the current to drop through zero to a 
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FIG. 3: Temperature and static force dependence of the flex- 
oelectric current for the 100/im sample in figure [2] As the 
static force is increased, the domain walls become less mobile 
due to impingement and a recovery of the positive flexoelec- 
tric current is observed. Similar behaviour was observed for 
the original thicker crystals. Domains with the c-axes along 
xi and X3 are labeled + and — respectively. 



negative value (bottom curves in figure [3]). The fact that 
the change of sign in the current is only seen when do- 
main walls can move leads us to believe that the domain 
walls may be charged with a polarization of opposite sign 
to that of bulk flexoelectricity, so that their motion un- 
der the periodic stress produces a current of the opposite 
sign to the flexoelectric response. 

The possible polarization of domain walls in STO (a 
non-polar material) is an unexpected result. At present, 
we can think of three explanations for the observed do- 
main wall charge: (i) the local strain gradient at the 
walls polarizes them through flexoelectricity, (ii) order- 
parameter coupling between the ferroelastic distortion 
and the (suppressed) ferroelectric polarization induces 
polarization at the domain wall [19| . or (iii) the domain 
walls trap charged deflects such as oxygen vacancies [2(| ■ 
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TABLE I: Orientation dependence of the flexoelectric re- 
sponse. 

We have also measured the flexoelectric response of 
[001], [101] and [lll]-oriented samples (inset of figure [2]) . 
The corresponding in-plane orientations were determined 
by X-ray diffration and are shown in table [H together 
with the average measured flexoelectric coefficients. All 
the coefficients were found to be of the same order of 
magnitude (ranging between about 1 and lOnC/m), but 
have different signs. While there was some inter-sample 



variation in the magnitudes of fj, for each of the orienta- 
tions, the signs were robust. 

For any material belonging to one of the cubic point 
groups there are only five independent components of 
the flexoelectric tensor fijki- In the case of STO, which 
belongs to the Oh group, the 4-fold rotation symme- 
try further reduces this number to three: fun, /1122 
and /i22i(= /i2i2)- For different crystal orientations, 
the measured polarization arises from different combina- 
tions of the three flexoelectric tensor components, i.e., 
the calculated values of n are effective coefficients rather 
than the flexoelectric tensor components defined in |2]). 
In addition we must not forget the contributions to 
from the gradients of e 2 2 and £33. For a bent plate [16| 
£33 = ~ §^ £ ii — f^f 622 wnere c ij are the elastic moduli. 
In our bending geometry, the anticlastic strain £22 is as- 
sumed to be negligible, giving for [001] oriented samples 
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The relevant anisotropic Poisson ratios z/y can be ob- 
tained from the known elastic moduli of STO F° r 
samples whose edges are not aligned with the crystal- 
lographic axes the fijki above should be replced by 

fijki = -LipLjqLkrLisfpqrs with Lij = X,; • &j and Vij by 

the corresponding v'^ . 

Inserting the values for /j, and v'^ for the different orien- 
tations of the crystals and the corresponding expressions 
for f-j k i in terms of fijki into equation ((4]) leads to three 
simultaneous equations. However, it turns out that these 
are not independent and hence cannot be solved to find 
the individual tensor components, leaving us instead with 
the following relations: 
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Nevertheless, the 



^ and 8 = C11 + C12 I 2C44 . 
interdependence of the equations can be used to check 
the self consistency of our analysis since we can use 
the measured /i values for the [001] and [101] to cal- 
culate the expected \x for [111] which turns out to be 
— 3.1nC/m and compares reasonably well, considering 
the intersample variation, with the actual experimen- 
tal value of — 2.4nC/m. A third independent equation 
can be obtained by changing the geometry of the sam- 
ple from a plate to a beam. In this case £22 = — Vi\t\\ 
and £33 = — ^3i£n [H}. Using this, rough estimates of 
/1111 ~ -9nC/m, /1122 ~ 4nC/m and / 12 i 2 « 3nC/m 
were obtained [22| . 

Taking these as order of magnitude values, we can es- 
timate the expected flexoelectric polarization around de- 
fects such as dislocations in epitaxial STO. Chu et al. 
and Nagarajan et al. report strains at dislocations 
of order 0.05 relaxing over several nm thus giving rise 
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to strain gradients of order 10 7 m _1 . Flexoelectric coef- 
ficients of order 10 _9 -10 _8 C/m will therefore give local 
polarizations of about 1-10/zC/cm 2 at room temperature, 
and up to 100 times more at low temperatures due to the 
increase in dielectric permittivity. This implies that the 
local polarization around defects in non-polar materials 
such as STO can be bigger than the ferroelectric polar- 
ization in the best ferroelectrics. Obviously, this is only a 
local effect, but given high enough density of dislocations 
(as can happen in strained thin films) we should expect 
the impact of flexoelectric polarization on the functional 
properties of dielectrics to be very large. 

Before concluding we should briefly address the possi- 
bility of artefacts. For high purity STO we do not expect 
any bulk piezoelectric contribution to the measured cur- 
rent. The lowering of the symmetry at the surface how- 
ever, introduces the possibility of a contribution from 
surface piezoelectricity @. Experimental data and first 
principles calculations suggest that the perturbed layer 
in only a few lattice constants thick, with no consensus 
as to whether it is polar or not 2J, |25j . Our estimates 
show that to mimic the measured response the piezo- 
electric modulus of the surface has to be approximately 
the same as that of a good ferroelectric such as BaTi03 
which seems unlikely. Finally, whatever the nature of the 
surface, there will in general always be a surface flexo- 
electric contribution f s f ~ e/a Q. However unlike the 
bulk effect, f s f is not expected to scale with the dielec- 
tric permittivity @, @] and thus it should be about two 
orders of magnitude lower than the bulk effect in STO. 
Nonetheless, recent works have revealed extended near 
surface skin regions up to 100/im deep with local fluc- 
tuations of the ferroelastic phase transition [26j| , Such 
regions were found to be inhomogeneously strained and 
thus may possess flexoelectric polarization even in the ab- 
sence of external stress. This may play some role in our 
samples, but even more so in fine-grained ceramics due 
to their higher density of surfaces, which may contribute 
to the very high values of the flexoelectric coefficient ob- 
tained in ferroelectric ceramics. Q- 

To conclude, we have measured the dielectric polar- 
ization induced by bending in single crystals of SrTiOa. 
Measurements of samples with different crystallographic 
orientations have allowed all components of the flexo- 
electric tensor to be estimated. These are of the order 
of 10~ 9 -10 -8 C/m, producing, around dislocations or de- 
fects, local polarizations of several /zC/cm 2 and higher 
at low temperatures. The analysis of the behaviour of 
the flexoelectric current as a function of static bending 
in the low-temperature phase also suggests that the do- 
main walls of STO are polarized, either intrinsically due 
to local gradient coupling, or extrinsically through defect 
accumulation. 
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